Edge states and topological phases in one-dimensional optical superlattices 
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We show that one-dimensional quasi-periodic optical lattice systems can exhibit edge states and 
topological phases which are generally believed to appear in two-dimensional systems. When the 
Fermi energy lies in gaps, the Fermi system on the optical superlattice is a topological insulator char- 
acterized by a nonzero topological invariant. The topological nature can be revealed by observing 
the density profile of a trapped fermion system, which displays plateaus with their positions uniquely 
determined by the ration of wavelengths of the bichromatic optical lattice. The butterfly-like spec- 
trum of the superlattice system can be also determined from the finite-temperature density profiles 
of the trapped fermion system. This finding opens an alternative avenue to study the topological 
phases and Hofstadter-like spectrum in one-dimensional optical lattices. 
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Introduction.- In recent years ultracold atomic systems 
have proven to be powerful quantum simulators for inves- 
tigating various interesting physical problems including 
many-body physics and topological insulators In 
comparison with traditional condensed matter systems, 
cold-atom systems provide more control in construct- 
ing specific optical lattice Hamiltonians by allowing both 
tunable hopping and trapping potential to be adjusted as 
needed, thus optical lattices populated with cold atoms 
offer a very promising alternative avenue to build topo- 
logical insulating states [2]. Several schemes exploring 
topological insulators with or without Landau levels in 
optical lattices have been proposed . So far all these 
schemes focus on two-dimensional (2D) systems, as one- 
dimensional (ID) systems without additional symmetries 
are generally thought as lack of topological nontrivial 
phases [?]. 

In this work, we study properties of trapped fermions 
on ID quasi-periodic optical lattices and show that these 
systems display non-trivial topological properties, which 
share the same physical origins of topological phases of 
2D quantum Hall effects on periodic lattices [8|. The 
quasi-periodic optical lattices can be generated by su- 
perimposing two ID optical lattices with commensurate 
or incommensurate wavelengths which has led to 

experimental observation of Anderson localization of a 
noninteracting BEC of '^^K atoms in ID incommensu- 
rate optical lattices [13] ■ Motivated by the experimen- 
tal progress, one-dimensional optical superlattices have 



been theoretically studied [I2| - |l6l |. For a quasi-periodic 



superlattice, the single-particle spectrum is organized in 
bands. When fermions is loaded in the superlattice, the 
system shall form insulators if the chemical potential 
(Fermi energy) lies in the gaps. For the open boundary 
system, localized edge states are found to appear in the 
gap regimes. The appearance of edge states is a signature 
indicating that the bulk states are topological insulators 
characterized by nonzero Chern number. We show that 



the topological invariant Chern number can be detected 
from the density distribution of a trapped fermion sys- 
tem, which displays plateaus in the local average den- 
sity profiles with positions of plateaus uniquely deter- 
mined by the ration of wavelengths of the bichromatic 
optical lattice. Through the analysis of universal scaling 
behaviors in quantum critical regimes of conductor-to- 
insulator transitions, we further display that both the 
positions and widths of plateaus can be read out from 
the finite-temperature density distributions, which pro- 
vides us an alternative way to study topological phases 
and Hofstadter-like spectrum by using ID optical super- 
lattices. We notice that localized boundary states in ID 
incommensurate lattices have been experimentally ob- 
served very recently by using photonic quasi-crystals ■ 
Quasi-periodic lattices.- We consider a ID polarized 
Fermi gas loaded in a bichromatic optical lattice [13, [HI , 
which is described by 

L 

= -t^(ctc,+i+H.c.)+^l^,n, (1) 

i i—1 

with 

V^ = Vcos{2Trai + S), (2) 

where L is the number of the lattice sites, cj (ci) is 
the creation (annihilation) operator of the fermion, and 
hi — c\ci. The hopping amplitude t is set to be the 
unit of the energy (t = 1), and V is the strength of 
commensurate (incommensurate) potential with a being 
a rational (irrational) number and 5 an arbitrary phase 
whose effect shall be illustrated lately. Suppose that the 
n-th eigen-state of a single particle in the ID lattice is 
given by \ipn) = '^iUi,nc\\0), the eigenvalue equation 
H\ipn) = En\ipn) Icads to the following Harper equation: 



(wi+i.n + Ui-i.n) + Vcos(27rQ;i -I- S)u, 



EnU, 



(3) 



where iti.„ is the amplitude of the particle wave function 
at i-th site with Vi the on-site diagonal potential and £"„ 
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FIG. 1: (Color online) Spectrum of Hamiltonian (1) with 
a = 1/3, V — 1.5 and 5 = 2n/3. (a) Energy bands for 
the system with 99 sites under periodic boundary conditions. 
(b)Eigen-energies in ascending order for the system with 98 
sites under open boundary conditions, (c) From up to down, 
figures represent the wavefunctions of states marked by labels 
("+" for states in bands and "*" for states in gaps) in (b). 



the n-th single particle eigenenergy. The ground state 
wave function of the TV spinless free fermionic system 
can be written as j^Pj?) = Yln=i Wi,nc||0), where N 

is the number of fermions. 

The solution to Eq. Q is closely related to the struc- 
ture of the potential Vi. For the incommensurate case 
Vi = Fcos(27rai) with irrational a, Eq.(l3]) is the well- 
known Aubry- Andre model (3], which has been showed 
that when y < 2 all the single particle states are ex- 
tended and when V > 2 all the single particle states are 
the localized states. However, for a periodic Vi, all the 
single particle states are extended band states according 
to Bloch's theorem. Next we shall consider the commen- 
surate potential Vi with a rational a given hy a ^ p/q 
with p and q being integers which are prime to each 
other. Since the potential Vi is periodic with a period 
q, the wave functions take the Bloch form, which fulfills 



Ui, for the lattice under the periodic bound- 
ary condition. Taking uj = e^^^cjijik) for |fc| < 7r/q, we 



have : 



(fc) 



In terms of 4>j{k), Eq.([3]) becomes 



[e'^(l)j+i + e ''"(jij-i] + Vcos{2njp/q + S)(f>j 



(4) 



Since ipj+q — 0j , the problem of solving the Harper equa- 
tion Eq.dU reduces to solving the eigenvalue equation, 
M$ = E^, where $ =(</>!, ••• , (fiqY and M \s a. q x q 
matrix. Solving the eigenvalue problem, we get q eigen- 
values: Ea{k) with a = 1, • • • Consequently, the en- 
ergy spectrum consists of q bands. As an example, the 
energy spectrum for a commensurate lattice with a — l/i 
and L = 99 under periodic boundary condition is given 
in Fig. la. For the lattice with open boundary condition. 



the momentum k is no longer a good quantum number. 
There appear edge states in gap regimes as indicated by 
symbols of "star" located at gap regimes in Fig. lb. As 
shown in Fig.lc, these edge states are localized in the 
left and right boundaries in contrast to extended states 
marked by symbols of in band regimes. 

As the phase 5 varies from to 27r, the spectrum for 
a given a = p/q changes periodically. The position of 
the edge states in the gaps also varies continuously with 
the change of 5. In Fig. 2 we show the spectrum of the 
quasi-periodic systems with a = 1/3 and 1/4 versus 5 
under the open boundary condition. The shade regimes 
correspond to the band regimes and the lines between 
bands are the spectra of edge states. 

Topological invariant.- Appearance of edge states is 
generally attributed to the nontrivial topological prop- 
erties of bulk systems 0, [3l ■ Next we explore the topo- 
logical properties of states under the periodic bound- 
ary condition. The topological property of the system 
can be understood in terms of Berry phase in k space, 
which is defined as 7 = § A^dk, where is the Berry 
connection A^ — i{(t>{k)\dk\4'{k)) and (/)(fc) the occupied 
Bloch state [20]. Adiabatically varying the phase 5 from 
to 27r, we get a manifold of Hamiltonian H{5) in the 
space of parameter 5. Similarly, we can define the Berry 
connection As = i{4){k,5)\ds\(t){k,6)) . For eigenstates 
<j}{k,5) of H{5), we may use the Chern number to char- 
acterize their topological properties. The Chern num- 
ber is a topological invariant which can be calculated via 
C — Jq^^"^ dk J^^ dS[dkAs — dsAk]. To calculate it, we 
need work on a discretized Brillouin zone. Here we follow 
the method in Ref. [2l| to directly perform the lattice 
computation of the Chern number. For the system with 
a = 1/3, we find that the Chern number for fermions 
in the lowest filled band (with 1/3 filling) is 1, while the 
Chern number for states with the second band fully filled 
(with 2/3 filling) is -1. 

To understand the topological origin of fermions in 
the ID quasi-periodic lattices, we explore connection of 
the present model to the well-known Hofstadter prob- 
lem [l^l , which describes electrons hopping on a 2D 
square lattice in a perpendicular magnetic field, with 
the Hamiltonian given hy H = ^ j) ^ii^]^*^*^'^'^'^ ' 
where the summation is over nearest-neighbor sites and 
the magnetic flux through each plaquette given hy (f) — 
Y^piaquette^ii- Taking the Landau gauge, the eigen- 
value problem is described by the Harper equation [8|: 
-tx{ijjj-i + V'i+i) - 2tyCOs{2iTj(j) - ky)iljj = E{ky)^j, 
where tx {ty) is the hopping amplitude along the x (y) 
direction. If we make substitutions oit^tx, V ^2tj,, 
a ~^ (j) and 5 — t- ~ky, the current ID problem can be 
mapped to the lattice version of the 2D integer Hall ef- 
fect problem. For the latter case, it was shown that when 
the chemical potential lies in gap regimes, the Hall con- 
ductance of the system is quantized [23*1 and the cor- 
responding Hofstadter insulating phase is a topological 
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FIG. 2: Energies varying with the phase 5 for systems with 
V = 1.5, (a) a = 1/3, L = 98, and (b) q = 1/4, L = 99 under 
open boundary conditions. 
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FIG. 3; The butterfly-like energy spectra with respect to a 
varying from to 1/2 with different phases: (a) 5 = 0, (b) 
6 = 7r/4. Both figures are for the system with V = 2 and 
L = 120 under periodic boundary conditions. 



insulator characterized by nonzero Chern number. 

Keeping this connection in mind, it is not strange to see 
that the energy spectrum of the ID quasi-periodic sys- 
tem versus different a has similar structure as the spec- 
trum of the 2D Hofstadter butterfly. In Fig[3^ and b, we 
show the spectrum for the system with phase 6 = and 
6 = 7r/4, respectively. The basic structure is quite similar 
for different phases except for some minor differences, for 
example, the spectrum for a = 7r/4 is continuous around 
i? = for (5 = 0, but separated for i5 = 7r/4 (see also 
Fig2.b). The familiar Hofstadter spectrum is actually a 
summation over all the ID spectra with different phases 
6 from to 27r. Despite the existence of the mapping 
between the ID quasi-periodic system and the 2D Hof- 
stadter system, we note that the edge modes in the ID 
system are localized and would not be used for dissipa- 
tionless transport as in high-dimensional systems. 



Experimental detection.- In realistic ultracold atom ex- 
periments, we need consider the effect of an external con- 
fining potential, i.e., Vi in Eq.([T|) is given by 



Vi = Vcosia2Tri + S) + Vh {i - io)^ 



(5) 



where Vh is the strength of the additional harmonic trap 
with io being the position of trap center. The den- 
sity distribution for the trapped system can be calcu- 
lated via Ui = (^'fI^iI^f) with ^I^p the ground state 
wavefuction. In Fig. 4, we have shown the local average 
density distribution of fermions trapped in both com- 
mensurate and incommensurate optical lattices with a 
harmonic trap. Here, in order to reduce the oscilla- 
tions in density profiles induced by the modulation of 
potentials, we have defined the local average density 



2-1-1=- 



■M 



i+j/i2M -hi), where 2M -f 1 is the length 



to count the local average density with M <^ L [2J]. Af- 
ter counting on the locations of plateaus, we find that 
heights of plateaus are completely decided by a with val- 
ues n(a) = a, 1 — a, 2a, 1 — 2a, 4a, 1 — 4a, if the values 
are in the range of (0, 1). The plateaus with fii = a, 1 — a 
are the widest ones corresponding to the largest gap 
regimes in the butterfly spectrum, while hi = 2a, I — 2a 
correspond to the smaller gap regimes. The width of 
plateau is associated with the size of gap. 

One can also understand the appearance of plateaus 
under the local-density approximation (LDA), in which 
a local chemical potential is defined as = fJ. ~ Vnii), 
where Vff (i) is the harmonic trap potential and /i is de- 
termined by X^i ^) = ^- The LDA is applica- 
ble provided that the number of fermions is large, and 
the variation of the trap potential is slow. Within the 
LDA, the local chemical potential ^{i) decreases parabol- 
ically away from the center of the trap. When the local 
chemical potential lies in one of the gaps, there appear 
plateaus in the density profile. The discernible num- 
ber of plateaus is related to the size of the energy gaps. 
For instance, in Fig. 4, the plateau with n = 1 corre- 
sponds to the band insulator with completely filled band, 
which is topologically trivial with zero Chern number. 
For a = 1/3 the chemical potential passes through two 
gap regions which produces two plateaus with n = 1/3 
and n = 2/3, whereas for a = 1/5 the chemical po- 
tential passes through four gap regions which gives four 
plateaus with n = 1/5,2/5,3/5,4/5, respectively. The 
Chern number can be obtained from the density by us- 
ing the Streda formula [2^ , which is valid when the 
chemical potential lies in a gap. From the Streda formula 
C = ^^Qa^ : can get C = 1, —1 for fij = a, 1 — a and 
C = 2, -^2 for Hi = 2a, 1 - 2a. 

In principle, all gaps in the spectrum, including the 
nontrivial smaller gaps in the butterfly, can be deter- 
mined via observing the corresponding plateaus. Widths 
of plateaus are proportional to sizes of corresponding 
energy gaps. Experimentally it becomes increasingly 



4 



Mott plateau 




100 200 j 300 400 500 



100 200 j 300 400 500 



FIG. 4: The local average density profiles for systems with 
rational a (a) and irrational a (b). Inserts for both pictures 
are the corresponding density profiles. The system is with 
1000 sites, 600 free fermions, V = 1.5, 5 = 0, Vh = 3 x 10"^ 
Here we take M — q for the rational case of (a) and M = 20 
for the irrational case of (b). 
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FIG. 5: (a): Local average density profiles for systems with 
different temperatures. Insert: enlargement of the plateau 
area, (b): Scaled local average density profles vs at dif- 
ferent temperatures. Insert: the universal function around 
critical point fj,c2 for the corresponding system. The system 
is with 1000 sites, 600 free fermions , V = 1.5, a — 1/3, 5 = 0, 
and Vh = 3x 10"^ 



harder to observe these smaller gaps, as the corre- 
sponding plateaus are very narrow which needs more 
precisely experimental measurement of the density dis- 
tribution. For the realistic detection, we need con- 
sider the effect of temperature. The finite-temperature 
density distribution can be calculated via ni{T) = 

|E^liie-^"/'"^(*?|clc.|*F>: where Ns = Ll/{L - 



Ny.N\ is the number of states, En is the energy of eigen- 

Ns ^-E„/kBT 



tition function 26|. In Figl5^, we display the local aver- 
age density profiles at different temperatures. It is shown 
that plateaus already become invisible for T > 0.3t. De- 
spite the fact that the obvious plateaus are smeared out 
by temperature fluctuations, we demonstrate that the 
position and the width of zero-temperature plateaus can 
be uniquely mapped out from finite-temperature density 
distributions of the trapped optical lattice system, which 
fulfil some universal scaling relations near the zero- 
temperature phase transition point fi = fic- 



^4-1 J- 



(6) 



state 4*^, and Z 



y 



is the canonical par- 



where n = n(T, /i) represents the density distribution 
for fermions with temperature T and chemical potential 
/i, nr{ii,T) is the regular part of the density, n(x) is a 
universal function describing the singular part of density 
near criticality, d = 1 the dimensionality of the system, 
v the correlation length exponent, and z the dynami- 
cal exponent. It was shown that 2 = 2 and ly = 1/2 
for the metal-insulator transition of ID free fermions 
(27} . From Fig. 5b, we see that the different curves of 
[n{fx,T) — 2/3]/T'^/^ intersect at two points. From inter- 
secting points, we can determine /Xc^ and /icj, and thus 
the corresponding gap in the spectrum can be inferred 
from /ic2 ~ Mci • In terms of the scaled chemical poten- 
tial /i = (/i — /ic)/T, curves for different temperatures 
collapse into a single one as shown in the inset of Fig. 5b. 

In order to connect to the real experiment in cold 
atoms, we refer to the parameters in Ref. (28| in which 
'^'^K atoms were used in an optical lattice with spacing 
a — 413 nm, the deep of the potential Vq = hEp, and 
recoil energy Eji = 45.98fi. kHz, thus the hopping mag- 
nitude t = 0.066Eji = 3.035fi, kHz. Correspondingly, 
the temperature T in Fig. 5 in unit t/ks, where fc^ is 
the Boltzmann constant, is of the order of 10 uK, e.g., 
T = 0.1 corresponds to 2.3084 nK. 

Summary.- In summary, we explored edge states and 
topological nature of Fermi systems confined in ID op- 
tical superlattices. Our study reveals that the topolog- 
ical invariant can be detected from plateaus of density 
profiles of the trapped lattice systems. Our results also 
clarify the connection of ID superlattice systems to 2D 
Hofstadter systems and will be useful for studying topo- 
logical phases and observing the Hofstadter spectrum by 
using ID optical lattices. 
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